Second-order Optical Response from First Principles 
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We present a full formalism for the calculation of the linear and second-order optical response for semicon- 
ductors and insulators. The expressions for the optical susceptibilities are derived within perturbation theory. 
As a starting point a brief background of the single and many particle Hamiltonians and operators is provided. 
As an example we report calculations of the linear and nonlinear optical properties of the mono-layer InP/GaP 
(110) superlattice. The features in the linear optical spectra are identified to be coming from various band com- 
binations. The main features in the second-order optical spectra are analyzed in terms of resonances of peaks in 
linear optical spectra. With the help of the strain corrected effective-medium-model the interface selectivity of 
the second-order optical properties is highlighted. 

PACS numbers: 61.50A, 42.65K, 61.50K, 78.66H 
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I. INTRODUCTION 

A material interacting with the intense light of a laser beam 
responds in a "nonlinear fashion". Consequences of this are 
a number of peculiar phenomena, including the generation of 
optical frequencies that are initially absent. This effect al- 
lows the production of laser light at wave lengths normally 
unattainable by conventional laser techniques. So the appli- 
cation of nonlinear optics (NLO) range from basic research 
to spectroscopy, telecommunications and astronomy. Second 
harmonic generation (SHG), in particular, corresponds to the 
appearance of a frequency component in the laser beam that 
is exactly twice the input one. SHG has great potential as a 
characterization tool for materials, because of its sensitivity to 
symmetry. Today SHG is widely applied for studying the sur- 
faces and interfaces because it requires an inversion asymmet- 
ric material. For materials with bulk inversion symmetry SHG 
is only allowed at surfaces and interfaces. This makes SHG 
a powerful surface selective technique. SHG in conjunction 
with Kerr and Faraday rotation 1 is used as an excellent tool 
for studying magnetic surfaces. In case of embedded inter- 
faces this technique gains extra weight when an intense laser 
is used which is capable of penetrating deep into the material 
and no direct contact with the sample is needed. 

In the case of linear optical transitions an electron absorbs 
a photon from the incoming light and makes a transition to 
the next higher unoccupied allowed state. When this electron 
relaxes it emits a photon of frequency less than or equal to 
the frequency of the incident light (Fig. 1(a)). SHG on the 
other hand is a two photon process where this excited elec- 
tron absorbs another photon of same frequency and makes a 
transition to yet another allowed state at higher energy. This 
electron when falling back to its original state emits a photon 
of a frequency which is two times that of the incident light 
(Fig. 1(b)). This results in the frequency doubling in the out- 
put. 

In order to extend the use of NLO for understanding the 
properties of surfaces and interfaces and for extracting max- 
imal information from such measurements for non centero- 
symmetric materials, a more quantitative theoretical analysis 
is required. The calculation of the SHG susceptibility from 
first principles is an important but difficult task. The major 
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FIG. 1: Schematic representation of a (a) linear optical transition 
and (b) second harmonic generation. 



work in this direction for semiconductors can be found in 
Refs. I2l3l4l5l and for metals in Refs. 16171 8l9l and references 
there in. In this work, we present the formalism for calcu- 
lating the second order susceptibility x^ 2 )(2w, uj, uS) for non 
magnetic semiconductors and insulators, within the indepen- 
dent particle approximation. The expressions formulated are 
amenable for numerical calculations using any band structure 
method and are computationally more efficient 10 than the pre- 
viously presented expressions. 2 £££ 

As an example, we present the linear and nonlinear opti- 
cal properties of a InP/GaP superlattice (SL). Semiconduct- 
ing strained SLs are potential materials for applications in op- 
tical communications involving switching, amplification and 
signal processing. In particular III-V semiconductor hetero- 
structures and SLs have attracted a great deal of interest 
mainly due to the possibility of tailoring band gaps and band 
structures by variation of simple parameters like superlattice 
period, growth direction and substrate material. Much of the 
theoretical work done to understand the physical properties of 
SLs has been largely concerned with the understanding of the 
electronic band structure. For example, the effect of strain on 
the band gap, the band offset problem and the possibilities of 
engineering it as well as the interface energy and band struc- 
ture have been studied. 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 The major 
theoretical work in the direction of NLO properties was done 
by Ghahramani et alr22&&2& They employed the non-self 
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consistent linear-combination-of-Gaussian-orbitals (LCGO) 
method to calculate the band structures and optical properties 
of SLs. As an example, we present a fully self-consistent cal- 
culation of the nonlinear optical properties of the mono-layer 
InP/GaP (110) superlattice (SL). The structure in the linear 
optical spectra is identified from the band structure of the ma- 
terial. The features in the second-order optical response are 
further analyzed in terms of resonances of the peaks in the lin- 
ear optical spectra. With use of the strain corrected effective- 
medium-model (SCEMM), 10 and we identify the features in 
the optical spectra of InP/GaP coming from the SL formation. 

The paper is arranged in the following manner. In Section 
II we present the detailed formalism for calculating the SHG 
susceptibility. Section III deals with the linear and second- 
order optical response of an InP/GaP SL. 



II. FORMALISM 

The aim of the present work is to obtain the expressions for 
the optical response of a material. The effect of the electric 
field vector E(w) of the incoming light is to polarize the ma- 
terial. This polarization can be calculated using the following 
relation: 



with 



P>)=xlV-EV)+X^ c -EV)-E 



(2) 



(1) 



In this expression x^ 1 ' is the linear optical susceptibility and satisfying 



%( 2 ) is the second order optical susceptibility. The higher 
order terms can also be calculated but in the present work 
we are only interested up to the second-order optical re- 
sponse. The formulae for calculating the SHG susceptibil- 
ity, x (2Wj w j w )> have been presented before. 2 - 3,4,5 In the 
present work we provide the detailed formalism. Before de- 
riving the expressions for calculating the SHG susceptibility 
using perturbation theory, a brief background of the single and 
many particle Hamiltonians and operators is needed. This is 
presented in the following three sections. 



A. THE HAMILTONIANS 

The optical response is treated within the independent par- 
ticle approximation and the Hamiltonian in a.u. (i.e. Tl = m = 
e = 1) is written as 



(2) 



The subscript i labels the electrons in the crystal at position Xj . 
p is the momentum operator given by pi = — iV». V(x) is 
the effective periodic crystal potential. K(t) = A(f)/c where 
A(i) is the vector potential of the external applied field. The 
macroscopic electric field is given by E(f) = — A(t)/c. In 
the long wavelength limit the variation of this field over the 
distance of the lattice spacing is neglected. Eq. (|2j can be 
separated into a time independent (which may be implicitly 
time dependent) and explicitly time dependent parts as 



Ho = h m = \ E p? + y ( x *) w 



J7i(t) = -K(t)£ Pi (5) 



H 2 (t) = -NK 2 (t) (6) 

N is the total number of electrons in the volume fl of the 
crystal. In the long wavelength limit Hi only introduces a 
time dependent phase factor for the wave functions and hence 
can be neglected. H\ can be treated as a perturbation. The 
eigenstates of Hq are given by 



H ip n (k,x) = w n (k)V>n(k,x) 
^„(k,x) = r!- 1 / 2 U „(k,x)e lk - > 



(7) 



Next we consider the time dependent single particle Hamil- 
tonian H(t) = i(p - K(i)) 2 + V(x). The instantaneous 
eigenstates of this Hamiltonian are 



&,(k,x) = fi- 1 /V(k + K(t),x)e 4k 



£T(t)$ n (k,x) = w„(k + K(t))Vi n (k,x). 



(8) 



(9) 



-0„(k, x) is implicitly time dependent via K(t). If an or- 
thonormal set ?/>„ (k, x) satisfies Eq. (|9) at time t, then the 
same equation is satisfied at time t + dt if 



ih—i> n (Kx) = E V>m(k,x)/w(k,t)E(i) (10) 



H — Hq + Hi + H2 



(3) 



for each n. If w m (k + K) = w„(k + K) then 

E(t).V mn (M)=0. (ID 
As shown Appendix A the above conditions require: 

Mmn (k + K,t) = . V "7^ I (12) 
tu) mn {k + K) 

where 

V mn (k,t) = J ^(k, X ) e - 4K - x hzV](^(k,x)) 

x e lKx dx (13) 

B. THE OPERATORS IN SECOND QUANTIZATION 

We now introduce Fermionic raising and lowering oper- 
ators a+ and a n satisfying the anti-commutation relations 
{a+,a+} = {a n ,a m } = and {a+,a m } = 5 mn . Similar 
commutation relations are satisfied by 6+ and b n which are 
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Fermionic raising and lowering operators, implicitly depen- 
dent on time via the wave function. The Hamiltonians in the 
second quantized representation are 



H = y"^„(k)q+ k a„ k 

nk. 

H(t) = 5>„(k + K)6+ k 6„ k 



(14) 
(15) 



A unitary transformation operator can be used for going from 
the implicitly time dependent operators b n and 6+ to time in- 
dependent operators a n and a+ as 



Here 



a n = Ub n U + 
a+ = Ub+U+ 

u = Y,\s><s\ 

s 



(16) 



(17) 



such that 1 5 > are states given by Eq. 0, \S > are states 
given by Eq. (JSJl and the sum runs over all the states. The 
importance of such a transformation will become clear in the 
next section (see for example Eq 1261 and the discussion after 
it). Since our final aim is to calculate the linear and nonlinear 
optical response of the material another useful operator is the 
current density operator J which is given by 



i 

nmk 

J' = ^ a„ k «mkV m (k + K). 



(18) 



(19) 



As shown in the Appendix B the current density is related to 
the response of the materials via polarization as 



J = J 



= JA 



dP 

~dt 



(20) 



here Ja is the intraband current density given by Eq. I66land 
P, the effective polarization of the material given by Eq. Q. 
In second quantization P can be written as 



P = ^ X & «k & ™kAw(k, t) 

nmk 

P ' = O £ a «k a "ikMnm(k, t) (21) 



C. THE PERTURBATION APPROACH 

1. SCHRODINGER TO INTERACTION PICTURE 

With the above information now we are ready to study the 
dynamics of a many-particle system described by the density 



matrix p, which is specified by the following equation of mo- 
tion 

ip = [H,p] (22) 

Working with the transformed operators p 1 = UpU + the 
equation of motion becomes 

ip' = [H' + H' dl p'} (23) 

where 

H' d = -fiP'.E(i) (24) 
= - q^ k a m k/w(k,t)E(f) 



nmk 

h' = UHU+ 



(25) 



= y^n(k + K)a+ k q nk 



nk 



In this equation a+ and a n are time independent and hence 
the commutators of W at two different times vanish. U is the 
unitary operater given by the Eq. ^] 

Until now we have been working in the Schrodinger pic- 
ture. Now we change to the interaction picture, because the 
equation of motion Eq. J22I involving total Hamiltoninan sim- 
plifies to the form given by Eq. d28l . where only the pertur- 
bation term of the Hamiltonian Hd, is involved. This can be 
done by the use of the following relations 



W(t) 

dW{t) 
dt 



exp 



H'(t')dt' 



iW(t)H'{t) 



(26) 



p = Wp'W + 
5 m k = Wa m kW + 

»mk(l) = drake mV ; 

~ X] a+ k a mk ^„ m (k,i)E(£).(27) 



H d = WH' d W+ = 



nmk 



Here f m (k, t) = J_ w m (k + K.(t'))dt'. As shown in Ap- 
pendix C Eq. J22I now becomes 

ip = [H d ,p] (28) 
Integrating this equation we get 



P = Po + - 

i 



(if 



[H d (t'),PoW 



[H d (t'),[H d (t"),p ]}dt'dt"... (29) 



oo J — oo 
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2. EXPECTATION VALUES OF OPERATORS 

Our main interest lies in the single particle operators like 
= E a+ k a mk B„ m (k) = ^ &+ k &mk@nm(k) 

nmk nmk 

6' = ^ k flmkBnm(k) 
nmk 

S = Q^ k Qmk6„ m (k + K) 

nmk 

The expectation value of any single particle operator 9 can be 
calculated using the relation 

< 9 >= Tr(p6) = Tr(p'Q') = Tr(p o 0) (31) 

where is 

If*. 

e = e(<) + - / [e(t),F < i(t / )]* / 



(30) 



— oo 

t fit' 



[[e(t), H d (t% H d (t")]dt' dt'' + (32) 



oo J — oo 



1 

+ WJ-c 

Using the following properties of the trace 

Tr(p a+ k a mk ) = S nm f n (k) 
Tr (po[ a tk a ™k a pk' a gk']) = 8u,U'8nq$mpfnm(k) (33) 

and combining Eqs. ( I29> . < I3 1 1 and d32t the first term becomes 
< 9 >(o) - E 6 nn (k + K)/ n (k) (34) 



n k 



and the second term is 

< © >(i) = i 



nmk 

C 

x e 

J —oo 



^ fnm (k)0 n?Ti (k, t)c mn ( ' ^ 

ll/m " (k '*'V m n(k,i')-E(i')^' (35) 



The higher order terms can be written in a similar manner. 
Taking E(t) = J2p E(w /3 )e" i "' 3 *, using the following iden- 
tity 



e s ^L(t) = 
and substituting 



d 




e s(t)d_ 


\L(t)] 


dt 


I S(t). 


dt 


[s(t)\ 


S(t) 




- wpt) 




L(t) 


_ b 
r'mn 






S(t) 


= i(bJmn 


- up) 





(36) 



(37) 



in Eq. d35l we get 

< e > (1) = i f nm (k)e nm (^ty— {Kt) 

nmk 

(k,t) 

i{(jj mn {n) - up) 

e i(v mn (k.,t')-uifft') 



X < e 



(38) 



<5< 



(w mil («') - Wg) 



E c (t')*' ^ E b ( W/3 ) 



Now expanding E{t') in frequency components E(uS) and us- 
ing Eq. fl36b for integrating the right side of Eq. fl38i one can 
write < 9 >(i)=< 6 >(i,o) + < B >(i,i) + where 



<0>(i,o) = /™m(k)6 nm (k,t) 



x e 



nmk 



(Wmn — ^/3) 



E b ( W/3 ) 



(39) 



< © >(!,!)= * E /nm (k)6 nm (k,t) 



nmk 



and 



_6_ 



< e > 



/ i mn(kj 



Wmn — W/3 — W 7 

i ^ +u -' )t E b (LJp)E c (u 7 ) (40) 



(2,0) 



E 



i 



Wmn — W/3 — W 7 
6 ,,c 



/» ; e nm (k, %m ZM? 
x E b (^)E c (w 7 ) 



(41) 



D. LINEAR AND SECOND ORDER SUSCEPTIBILITY 

Up to this point the things have been most general. Now 
in order to find the linear and non linear susceptibility is 
replaced by the polarization operator P. Eq. ^ can be written 

as 



< P >=< P >/ + < P > 



ii 



(42) 



where 



<P a >// = Xabc '{-Up,-W 1 ,Up,U 1 ) 

x e- i ^ +w ^ t E b (wp)E c (u 1 ) (43) 

Similarly expanding the intraband current density in the pow- 
ers of the electric field we get 



with 



< J A >=< Ja >i + < Ja >ii 



<Ja>/ = <7$(-W0,W0)e-^»*EV/j) 
< Ja >ii = a^(-up,-u y ,wp,u;y) 

x e- i( ^+^ )t E b (w /3 )E c (w 7 ). (45) 

For clean semiconductors with filled bands the conductiv- 
ity c' 1 ' is zero and only Pj contributes to the linear term, 



(44) 
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whereas to the second-order term both, Pji and J//, con- 
tribute. Now substituting O = P a in Eq. ( 1391 we get the 
expression for the linear susceptibility 



-y 

nmk 



x$(-<»,<») = ^>;/ «** 



(k)r^n(k) 
,„(k) - u> 



£q&M - Sab 

4tt 



(46) 



Here e ab (uj) is the ab component of the dielectric tensor, and 
r nm are the position matrix elements given by 

V„ 



f^nm 

o 



,(M) 



,(k + K) 



if LO n ^ L0 n 



if Ul r , 



(47) 



Similarly, substituting = P a in Eqs. J40i and (14- 1 i and 
using the following identity 



7 (77) 

abc 



W/3 



" 7 '^'^ 8) 



7 

„3.10 



we get after some rearrangement of the terms • the fol- 
lowing contributions to the SHG susceptibility: the inter- 
band transitions Xi ntcr (2ti;, u>, lo), the intraband transitions 
Xintra(2w, ui, uj) and the modulation of interband terms by in- 
traband terms Xmod^w, u>, o->)'- 



Wk is the weight of the k point, n denotes the valence 
states, m the conduction states and / denotes all states (I ^ 
m, n). We have used the these expressions to calculate the 
total susceptibility in the following example. 



III. EXAMPLE 

As an example we present the linear and nonlinear optical 
spectra of an InP/GaP (110) superlattice. This material is a 
mono-layer SL in (110) direction with GaP grown on top of 
an InP substrate. The zz component of the linear frequency 
dependent dielectric function is given in Fig. 2(a). ej^o;) nas 
major peaks at 2.3eV (B), 4eV (D) and 5.5eV (F) and minor 
peaks at leV (A), 2.75eV (C), 4.5eV (E) and 6.5eV (G). These 
peaks in the linear optical spectra can be identified from the 
band structure. The calculated band structure along certain 
symmetry directions is given in Fig. 3. It can be noted from 
the band structure plot that InP/GaP is a direct band gap (Eq) 
material. The calculated band gap using the local density ap- 
proximation (LDA) is Eg = 0.6eV— . As can be seen from 
the Eq. J46> . in order to identify of these peaks we need to 
look at the optical matrix elements r nm for various pairs of 
band n and m. We mark the transitions, giving the major 
structure in e| z (w), in the band structure plot. These transi- 
tions are labeled according to the peak labels in Fig. 2(a). 



nm/k 



Total 

- - Inter (2d)) 

Inter (to) 

■ - ■ Intra (2d)) 
Intra (H>) 



X£» c ra (2u,,u,, w ) = 

x [ w i™rni{ r Zm r mn} ~~ u rnl*lm{ r mn Y nl}\ 
1 



8/ 




ran 



2cj) 



nm I ml In 1 



2 4 6 

Energy(eV) 



1 2 3 4 5 6 7 
Energy(eV) 



, T nm{ V mi r in}{ u> rnl ~ Wn) 

~, wa m (« mn -2w) 



(50) 



k 



U! mn (u! mn — U>) 
nml " LU ' 



x [w r i;r" m {r^„r^;} — w/ m r° ; {rf TO r^ ln }] 

a f b 1 ~| 

nm I mn mn J I 



Umni^r, 



W) 



(51) 



FIG. 2: (a) Imaginary part of the zz component of the linear di- 
electric tensor, (b) Second-order susceptibility Im [x « (2 Wj U , a,)] 
(solid line) and different contributions to it: the u 2u interband term" 
(dashed line), the "lj interband term" (dotted line), the "2oj intraband 
term" (dash dotted line) and the "w intraband term" (circle thin line) 

We now go on to study the NLO properties. Different con- 
tributions to the imaginary part of Xxyz{2ijJ,ui,ui) are pre- 
sented in Fig. 2(b). As can be seen the total SHG suscep- 
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(1/2,0,0) (1/2,1/2,0) (0,1/2,0) (0,0,0) (0,0,1/2) 



FIG. 3: Band structure for mono-layer InP/GaP SL. 
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FIG. 4: Full SL calculation (solid line) along with the SCEMM 
results (dashed line) for (a) the imaginary part of the linear dielectric 
tensor and (b) the magnitude of the SHG susceptibility 



cause of the complicated resonance of the 2lj and lu terms. But 
one can make use of the linear optical spectra to identify the 
different resonances leading to various features in the NLO 
spectra. This analysis is performed in the present work. The 
identified peaks are marked in Fig. 2, where the nomenclature 
adopted is M(xw) + N(yui), which indicates that the peak 
comes from an xui resonance of the peak M with the ycu res- 
onance of peak /V in the linear optical spectra. For example, 
the hump just below 1 eV, labeled A(2oj) in the imaginary 

(2) 

part of Xxyz(2uj,uj,uj) comes from the 2u> resonance of the 
peak labeled A in the linear optical spectra. The peak labeled 
D(2u>) + B(lo) is coming from the 2uj resonance of the peak 
D with uj resonance of the peak labeled B in the e2(uj) plot. 

Compared to the linear optical, the NLO is a much more 
surface/interface sensitive technique. This fact can be demon- 
strated by identifying the features coming from the inter- 
face formation. These features are referred to as SL fea- 
tures. These features can be pin pointed by comparing 
the spectra for the SL with features appearing from the av- 
erage of the two bulk materials. In order to provide a 
simple model for predicting the averaged bulk features in 
the optical properties on basis of its constituent materials, 
the effective-medium-model (EMM) and the strain-corrected- 
effective-medium-model (SCEMM) 10 have been proposed. 
Comparison of the SCEMM results with the SL calculations 
are presented in Fig. 4. The SL features coming from effects 
like symmetry lowering are not accounted for by the SCEMM 
and are marked as SLX in the figure, with X representing the 
feature label. The small SL effects in the linear optical spec- 
tra are greatly enhanced in the second-order optical response. 
This clearly indicates the selective interface sensitivity of the 
NLO. 
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IV. APPENDIX A 



To show Eq. [21we need to show: 
d 



dt 



[(H(t) - w„(k + K))^ n (k, x)] = (52) 



tibility is zero below half the band gap. The 2u> terms start 
contributing at energies ~ \E>g and the co terms for energy 
values above Eq. In the low energy regime (< ieV) the SHG 
optical spectra is dominated by the 2uj contributions. Beyond 
3eV the major contribution comes from the u> terms. Unlike 
the linear optical spectra, the features in the SHG suscepti- 
bility are very difficult to identify from the band structure be- 



which can be written as 

[H(t) - co„(k + K; 
d 



dtp n (k,x) 
1 dt 



dt 

eE(t) 



[#(*)- w„(k + K)]^(k,x) 
Su> n (k + K) 



6k 



(iV - K) 



(53) 
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substituting Eq. Eq. [53]can be written as 

- * ^ W „(k + K) /lmll .E(t)(i(k,x) e ,K ' x ) (54) 

m 

"(W„(k + K) 



further using the property of trace and operator relations 

Tr([A,B]C) = Tr(B[C,A}), (62) 



= eE(t). 



( 5k 



iV 



(^(k,x)e lK - x ) 



The prime indicates that m ^ n. Note that here we have used 
the fact iV(^„e iK ' x ) = e lK x (iV - K)V>„. The projection of 
LHS of this Eq. on ^(k, x)e iKx is 

- / </5;(k,x)e~ lK ' x i ^w m „(k + K)/w 

m 

. E(t)VAn(k,x)e iK - x dx (55) 

For the case with I = n (m ^ n), this is zero. The RHS of 
Eq. [33] is also zero because 

i y C(k,x) e - lK - x V^„(k,x)e lK - x dx = -V nn (k,t)(56) 

Su n (k + K) 
8k 

For the case / ^ n the LHS of Eq. I54l is 

- i J ^(k,x)e- 4K - x ^ WmnAi mn-E(i)(^m(k,x)e lK - x ) 

rn 

= - iw In (k + K) Wn .E(t) (57) 

The RHS of Eq. |54| in this case is -eVi n .E(t). Thus the 
condition for two sides to be equal is 

~^mn (k, t) 



H' = ^^(k + K)a+ k a mk , 



(63) 



nk 



[ a «k a "ik, H'} = -uj nm (k + K)a+ k a mk , (64) 



and 



(65) 



the first term in Eq. l61lbecomes 



p mn (k + K,t) 



iuj mn {k + K) 



(58) 



V. APPENDIX B 



In this section we confirm Eq|66| We start with the relation 

< P > = Tr(pP) = Tr(p'P') 



-Tr 



P' Y^ a n\s. a mk.P>nr 



(59) 



Now taking the time derivative of the expectation value of P 
we get 



d< P a > 
dt 



1 / d 
-Tr — 
ft \ dt 



P' Y a nk a mkPr: 



= -Tr 

n 



dp' \ - + 

1 1 / & n k ^ m k M n m 



nmk 



nmk . 



dt 



(60) 



Using the equation of motion[28]this can be written as: 

d < P a > 



dt 



\ nmk / 

+ ^}Tr ( [H' d , p'} J2 "ikWL 

\ nmk / 

n [ p ^ v ™ k ^ hr (61) 

\ nmk / 



Y a tu. a m^Pnm - Y a+ k a nk uJn J 

.nmk nk J / 

= i | -Tr L' £ V«„ (k + K)a+ k a„ k ^ 

+ Tr ( pi Y, < k «mkK<5„ m - MSm-^W) ) I (66) 

\ nmk / ) 

where the last term in this equation is Ja 
VI. APPENDIX C 

We need to prove 

iP=[H d ,p] (67) 
The LHS can be written as 

Vp = i±p = i±\W(/W+] 

= i^p' W + + lW d -fw+ + l Wp'°^ (68) 
dt dt dt 

Now using Eqs. |23]and[26]this becomes 

- WH'p'W + + W[H' + H' d , p']W + + Wp'H'W+ 

= W[H' + H' d ,p']W + -W[H' lP ']W + 

= W[H d ,p']W + = [H d ,p\. (69) 
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